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ABSTRACT 

The discovery that Vega is a rapidly rotating pole-on star has raised a number of questions 
about this fundamental standard, including such issues as its composition, and in turn its mass 
and age. We report here a reanalysis of Vega's composition. A full spectral synthesis based 
on the Roche model derived earlier from NPOI interferometry is used. We find the line shapes 
in Vega's spectrum to be more complex than just flat-bottomed, which have been previously 
reported; profiles range from slightly self-reversed to simple "V" shapes. A high SNR spectrum, 
obtained by stacking spectra from the ELODIE archive, shows excellent agreement with the 
calculations, provided we add about lOkms -1 of macroturbulence to the predicted spectra. From 
the abundance analysis, we find that Vega shows the peculiar abundance pattern of a A Bootis 
star as previously suggested. We investigate the effects of rotation on the deduced abundances 
and show that the dominant ionization states are only slightly affected compared to analyses 
using non-rotating models. We argue that the rapid rotation requires the star be fully mixed. 
The composition leads to masses and particularly ages that are quite different compared to what 
are usually assumed. 



Subject headings: line: profiles — stars: abundances 
— stars: individual (Vega) — stars: rotation 



stars: chemically peculiar — stars: early-type 



1. Introduction 

With the announcement of the detection of the 
i nterferometric s ignatu r e of rap i d rotation in Vega 



ap] 

( Peterson et all 12004 l2006bt lAufdenberg et al 



2006), a number of questions were raised about 
the fundamental standard. Earlier suggestions 
of rapid rotation wer e based on the high lu- 
minosity (|Petriel Il964 iGravl Il988h of the ob- 
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ject and the un usual shapes of the weak lines 
in the spectrum (iGulliver. Hill, fc Adelmanl 19941 
Hill. Gulliver, fc Adelmanl 12004) . The high lu- 



minosity is immediately explained using Roche 
models for the fi gure of the rota ting star, von 
Zeipel's theorem (jvon Zeipell 1 1924 ) to character- 
ize the temperature distribution, and adopting a 
nearly pole- on geomet ry, required by the small 
line widths (|Gravlll988f) . 

But this model, a star rotating near breakup, 
raises the question of whether such fundamen- 
tal issues as Vega's composition, mass, and 
age, are accurately known. It has been recog- 
nized fo^_some_Jhne_Jhat Vega appears metal 



poor ( Sadakane fc Nishimuralll981 : Adelman fc Gulliver 
1990). And although it has been known since the 
early 20th century that masses deduced from lu- 
minosity and radius measurements are strongly 
affected by composition, recent mass and age de- 
terminations have largely assumed solar compo- 
sition, the assumption being that sharp lined A 
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stars often show abundance peculiarities that are 
assumed due to diffusion and generally confined to 
surface layers. The recognition of rotation veloc- 
ities approaching breakup renders that assump- 
tion unlikely, since rotation driven circulation is 
likely to mix the envelope completely and deeply 
over times short compared to operable diffusion 
timescales. 

Furthermore, the large surface temperature 
gradients that would be associated with high ro- 
tation raise a new question: how seriously are 
simple, single model atmosphere analyses of the 
spectrum affected by the composite nature of the 
atmosphere? The peculiar line shapes add to this 
concern. A full analysis of the spectrum, or at 
least representative spectral features, seems nec- 
essary to demonstrate that we understand the pe- 
culiar line profiles and are able to derive reliable 
abundances. 

This in turn requires high resolution, low 



1994 iHill. Gulliver, fe Adelmanll2004f ). These au- 



noise; spectra comparab l e to 



Hill. Gulliver, fc Adelman (2004 



those used by 
The spectra 
of Vega available on the ELODIE archive provide 
us the necessary resolution and low noise, we de- 
scribe those data in §2. The computation of the 
synthetic spectra based on a rotational model are 
described in § 3 and the deduced abundances and 
other characteristics we reported in § 4, includ- 
ing the discovery that significant macroturbulence 
must be adopted. In § 5 we discuss the implica- 
tions of the abundance profile and argue that the 
suggestion that Vega belongs to the A Bootis class 
of objects is probably correct. We note that the 
effect of rotation on the line strengths depends 
strongly on the line considered and propose a sim- 
ple resolution to the prediction of large departures 
from LTE in the Fe i spectrum that have not been 
seen in practice. We examine more closely the is- 
sue of rotational mixing and conclude that the 
abundances we find here likely represent the ma- 
terial out of which Vega was formed. Lastly, we 
estimate the mass and age of Vega based on this 
composition. 

After submitt ing this manuscript we became 
awar e of a paper (jTakeda. Kawanomoto. fc Ohishi 



2008) that had been recently accepted for publica- 



thors draw a number of conclusions in agreement 
with what we find here. But they also arrive at 
quite a different physical model of Vega, conclud- 
ing in the process that errors were made in the 
reductions of the published interferometry. We 
will comment on these results at the appropriate 
points. 

2. The Observational Data 

The V ega spectra we used are from the ELODIE 
archive ( Moult aka et al.l l2004h . which contains 
high-resolution (R ~ 42, 000) echelle spectra from 
the ELODIE spectrograph obtained at the Obser- 
vatoire de Haute-Provence 1.93 m telescope. The 
ELODIE data pipeline automatically extracts the 
spectra, establishes the dispersion, and corrects 
for scattered light. The spectra used here were 
obtained between 1996 and 2004. The wavelength 
rectified spectra covering AA 4000-6800 are pro- 
vided with 0.05 A sampling. 

Barycenter corrections were required before co- 
adding the spectra. To improve the signal-to-noise 
ratio (SNR) which is 250 for a typical spectrum, 
we co-added 49 out of 71 available spectra. In 
the process we rejected spectra whose SNR was 
less than 100 and those showing noticeable fring- 
ing. We also replaced bad pixels whose residu- 
als in individual spectra were 5 times larger than 
the typical noise by interpolating adjacent pixels. 
The co-added spectra were converted to a residual 
intensity scale by normalizing them to the scale 
of the synthetic spectra described in §3. The re- 
sultant spectra, segments of which are shown in 
Figures Q] and [2J were then compared to the syn- 
thesized spectra for the abundance analysis. The 
SNRs of the co-added spectra were estimated to 
range from 750 to 2,200 depending on the spectral 
regions. 

Besides the ELODIE spectra, spectra of com- 
parably high SNR and resolution of Vega have 
been obtained at the Dominion Astrophysical Ob- 
servatory (DAO; iGulliver Hill fe Adelmanl Il994 ; 

IHill. Gulliver, fc Adelmanll2004l ) and the Okayama 

Astro physical Observatory (OAO; lTakeda. Kawanomoto. fc Ohishi 
l2007t ). The DAO spectra (SNR ~ 3, 300) have not 



tion in this Journal which u ndertakes an analysis , 

. , , , j /Ittt - i — i — Tr — o yyu been released publicly. However tficy are avail- 

of the ir previously published ( lakeda, Kawan omoto. fc Ohisnil J 



20071 ) spectra toward understanding Vega's rota- 



tion, much along the lines taken earlier (jGulliver. Hill, fc Adelman 
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Fig. 1. — Shown are line profiles representative of the range of shapes encountered for weak lines in the 
ELODIE spectra of Vega. The shapes run from weakly "self-reversed" (e.g., Fe i A 4528 and Ba n A 4554) 
through flat-bottomed (Cr n A 4565 and S i A 6052) to "V" -shaped (O i A 6046). Where known, blends are 
indicated in parenthesis. Wavelengths are in the star's rest frame. 
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Fig. 2. — Plotted here are additional segments of spectra (dotted lines) showing the range of shapes of weak 
lines, as in Figure [1] only now overplotted with the synthetic spectra (continuous lines). Note particularly He 
I A 4713 which, with an excitation potential of 21 eV, is formed in a small region around the rotational pole 
and displays the corresponding "V" shape. At the other extreme Ca i A 6162 shows the weak double-horned 
( "self- reversed" ) shape reflecting its very low excitation potential, 1.9 eV; it is contributed exclusively by the 
cooler equatorial regions. Other lines showing this behavior are Ti n A 4708, Fe i A 5586, and Ca i A 5588, 
although all three are (weakly) blended. Two iron lines, Fe n AA6147 and 6149, at intermediate excitations 
of 3.9 eV above the 7.8 eV ionization potential of Fe i, show the expected flat-bottomed shapes, although 
seen against a slight variability in the background continuum. The weak Ca i lines indicated with ":" were 
not included in the abundance determination. 
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able as part of a graphic toolkit^ which allows one 
to examine sections of spectra at high resolution 
and identify lines and probable blends. We made 
extensive use of this tool during this investiga- 
tioni 

Also recently published are OAO spectra (SNR 
from 1,000 to 2,000 on average) covering AA3900- 
8800. However, these spectra display emission (e. 
g.,~ A 4560) and absorption (e. g.,~A6060) fea- 
tures and show the head of the Paschen continuum 
to be strongly in emission, features not reported 
elsewhere. So we have chosen to focus exclusively 
on the ELODIE data set. 

3. Computations 

We assume Vega can be described by a gravity- 
darkened Roche spheroid in solid-body rotation, 
with a point mass gravitational potential, show- 
ing a temperature distributio n varying accord 
ing to von Zeipel's theorem ( von Zeipel 1924) 



and s e en nearly pole-on (e. g., iPeterson et al 



2006bl : lAufdenberg et al.ll2006l ). Because the re- 
cent interferomctric measurements taken at the 
Navy Prototype Opti cal Interferometer (NPOI; 



Armstrong et alj I1998T ) and th e Center for High 



Angu lar Resolution Astronomy (jten Brummelaar et al 
2005) array yield closely similar model parame- 
ters, we ado pt the parameters obt ained from the 
NPOI data ([Peterson et~aD l2006bh for synthesiz- 
ing spectra; the model has a fractional rotation 
velocity, u> = 0.926, a polar surface gravity of 
logg p = 4.074, a polar effective temperature of 
T p = 9988 K, an inclination of the rotational axis 
to the line of sight i = 4.54°, and a projected 
rotational velocity of wsini = 21.7kms _1 . For 
det ails of fitting Roche m odels to the NPOI data 
see 



Peterson et al.1 (l2006aD . and fo r issues specific 



to Vega see Peterson et al. ( 2006bl ). respectively. 



To calculate the emergent spectrum we con- 
structed a square 256 x 256 grid which contains the 
apparent disk of the star, calculated the stellar pa- 
rameters at the center of each cell that actually fell 
on the flattened disk, and computed an emergent 



J http: //www. brandonu. ca/physics/gulliver/ccd_atlases .html 

2 At the same time R.L. Kurucz (2007, private communica- 
tion) provided a high resolution synthesized spectrum for 
Vega in the 450—500 nm region based on a line list and 
gf values calibrated to a solar spectrum which proved ex- 
tremely useful. 



flux as a function of A, /i (cosine of the angle be- 
tween the local normal and the line of sight), T e ff, 
<7eff (local gravity reduced by centrifugal force), 
and projected ve locity using the ATLAS 9 model 
atmosphere grid ( Castelli fe Kuruc3l2003l ) and the 



at omic line data g i ven in the extensive compilation 
of iKurucz fe Belli (|l995l ). The fluxes were inte- 



grated over the disk to yield the synthetic model 
spectrum. In these calculations LTE, hydrostatic 
equilibrium, and plane-parallel atmospheres were 
assumed to represent the star's surface locally. 



A concerns have been ra ised recently (|Aufdenberg et al 

120061 : iMonnier et al. 2007 ) about the rigorous ap- 
plicability of the von Zeipel theorem in the parts 
of the disk of a rotating star that are rendered 
cool enough to generate convection. We believe 
the issue is not relevant to Vega. In our model the 
temperature drops to about 7600 K at the equa- 
tor, and the effective gravity in turn decreases to 
about logg ^3.5. From a model atmosphere with 
T c ff = 7500 K, log<? — 3.5 we find the reduced 
density and in turn increased fraction of hydrogen 
ionized compared to the main sequence, substan- 
tially decreases the extent of the convective region 
and the efficiency of the resulting convection. Con- 
vection carries significant flux only in the range 
of 1 < TR OSSe i an d < 30, well out from the interior 
where the flux requirement is established. 

4. Results 



4.1. Line Shapes 

The abundance analysis was done by adjusting 
each element abundance until the model spec- 
tra fit the co-added spectra. Since Vega's lines 
are sharp and blending is minimal, the pro- 
cess of adjusting the abundances was straight- 
forward. Several representative regions of the co- 
added ELODIE spectrum are shown in Figures 
Q] and [2] Weak lines throughout the spectrum 
show not only the flat-bottomed shapes (Cr 11 
A 4565, S 1 A 60 52, and Fe 11 A 6147) as noted hi 
recent studies (iGulliver. Hill, fc Adelmanl [l994; 



Hill. Gulliver, fc Adelmanl l2004h but also weakly 



"self-reversed" shapes such as Mg 1 A 4702 and Ca 
1 A 6162 and "V" shapes such as Hei A 4713 and 
O 1 A 6046. 

The unusual shapes of the weak lines are 
strongly correlated with excitation and ioniza- 
tion potential and can be understood in terms of 



5 



how the Boltzmann factors amplify the tempera- 
ture gradient across the disk. Since Vega is seen 
nearly pole on, the center of the apparent disk is 
almost exactly at one pole, the hottest point on 
the star. On the other hand, the limb is nearly 
the equator which is not only 2,400 K cooler than 
the pole, but the visible gas is actually cooler still 
owing to the simple projection effects associated 
with limb-darkening. Therefore the bound states 
responsible for the lines seen from the light ele- 
ments such as He i, O i, Mg n, Al n, and Si n whose 
ionization and excitation potentials are quite high 
are excited mostly at the axis with zero projected 
velocity. There is almost no contribution to the 
line profiles from the rotationally shifted equato- 
rial region, resulting in "V" shapes. The lower 
the excitation potentials the lines have, the more 
enhanced the contribution from the more rapidly 
rotating equatorial regions becomes and the wider 
and more square shaped the line profiles get. For 
the elements such as Ca i, Fe i, and Ba n with the 
lowest excitation potentials one sees a mild double- 
horned shape ( "self-reversed" ) as the contribution 
from the equatorial region completely dominates 
the profile. In this sequence the flat-bottomed 
shape is formed at intermediate excitation po- 
tentials such as those of the lines of Cr n and 
Fe ii. Our synthetic spectra predict well this se- 
quence of line shapes as shown in Figure [2j where 
three regions of the ELODIE spectrum, overplot- 
ted with our synthetic spectrum, are shown. We 
see that weak Fe n lines tend to have a flat- 
bottomed shape while weak Fe i lines show a self- 
re versed shape. Examinatio n of d ata presented 
bv lGulliver. Hill, fc Adelmanl lfl99l and particu- 
larly in Iffill. Gulliver, fe Adelmanl (j2004h suggests 
these shapes are present in their data as well. 

4.2. Macroturbulence 

In the process of the spectral synthesis, in or- 
der to fit the shapes of the weak line we found 
we had to reduce the resolution of the spectra 
well below the nominal resolution of 42,000 of the 
ELODIE spectra, ultimately adopting a resolution 
of about 25,000 as shown in Figure [31 We inter- 
pret this additional broadening, which was accom- 
plished by convolving the synthetic spectrum with 
a Gaussian, as adding lOkrns -1 of macroturbu- 
lence to the nominal ELODIE resolution (also as- 
sumed to be a Gaussian). The effect of this addi- 



tional broadening is most noticeable in steep-sided 
line profiles (e. g., Mg i A 4703 and Ni i A 4713), 
as shown in Figure [3l As a result the value de- 
termined for the macroturbulence comes from low 
excitation lines and hence refers more to the equa- 
torial regions than the polar regions. As might be 
expected, "V" shaped lines such as He i A 4713 are 
insensitive to the macroturbulence, as also shown 
in Figure [3l This is a very interesting result which 
we discuss at more length below. One caution is 
immediately apparent though: line widths might 
not be reliable indicators of actual projected ve- 
locity, at least for stars seen at low inclination. 

At this point our analysis deviates sharply from 



the re cent contribution from lTakeda. Kawanomoto. fe Ohishi 
(|2008l ). who seem not to have considered the pos- 
sibility of large scale non-thermal line broadening. 
That there could, and even should, be turbu- 
lence on large scales in the atmosphere of Vega 
seems easy to justify. Even very slow, cms -1 , 
subsurface circulation currents will be magnified 
by the many order of magnitude drop in den- 
sity found in the outer envelope, as required by 
the equation of continuity. Add to this a very 
strong Coriolis force owing to the rapid rota- 
tion and a surface covered with large eddies - 
cyclones - is to be expected. Ignoring t hi s pos - 
sibility, iTakeda. Kawanomoto. fe Ohishil (|2008h 
were forced to adopt a relatively slowly rotating 
model, creating a clear c onflict with the inter- 
ferometric measuremen ts (Pete rson et al . 200Gb; 
Aufdenberg et al.ll2006l) . 



4.3. Abundance Analysis and Microturbu- 
lence 

As is often the case, we found that it was 
generally not possible to find abundances for ele- 
ments (or even the same ion of an element) which 
gave good fits to both strong and weak lines si- 
multaneously. This is usually taken as a sig- 
nal that some microturbulence needs to be intro- 
duced. To this end, we determined the abun- 
dances from the Fe n lines for two choices of 
the microturbulence, as shown in Figure [U Here 
the abundances are given as the logarithm of 
the ratio of the number of an element to that 
of total elements 



Castelli & Faraggiana 
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(|1979l ) and ISadakane fe Nishimural (|198lh have 
previously noted that the influence of the micro- 
turbulence is less important in the visual region for 
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Fig. 3. — These figures plotted as in Figure [2] illustrate the need for line broadening in addition to rotation 
and microturbulence (in both panels the lower spectra are offset by 0.1). Mg i A 4702 shows the problem 
most clearly although it is also evident in Ni i A 4714. The nominal ELODIE resolution of 42,000 (assumed 
to be Gaussian) allows too much structure in the steep-sided line profiles. Reducing the resolution to 25,000 
appears to be required, which we interpret as a contribution of about 10 kins -1 of macroturbulence. The 
effect of adding this macroturbulence is to improve the fit dramatically in the bottoms of the weak, low- 
excitation lines while causing the line widths to be a bit wide. This suggests the actual projected rotation 
rat e is below the a d opted 2 1.7 kms -1 , as was suggested in the initial interferometric data reductions reported 
by IPeterson et al. (|2006tJl . Note that high excitation lines like He i A 4713 are not affected by the added 
macroturbulence. 
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lines of intermediate strengths about 40 ~ 70 mA , 
which we also found. For Fe n, which has the 
widest range of equivalent widths, we find both 
the scatter and any trend with equivalent width 
are significantly reduced for a microturbulence of 
about 2 kms -1 which we subsequently adopt. The 
O i triplets also support 2kms~ 1 (e. g., Figure [2] 
which shows only the case for 2kms~ 1 ). 

Table [1] shows the deduced abundances for Vega 
with a microturbulence of £p = 2kms _1 . In se- 
lecting lines we eliminated severe blends but in- 
cluded weak blends where we felt reliable abun- 
dances could be obtained. The columns are the 
laboratory wavelength, lower excitation potential, 
equivalent width, log gf, and the deduced abun- 
dance (log jj* 1 - ) . Blends we have decided to retain 
are noted in the last column. The abundances for 
elements with only single lines such as Al n, S i, 
Mn i, and Ni i must be considered uncertain. Even 
where there was no obvious blending, abundances 
were determined exclusively by spectral synthe- 
sis. Nevertheless, we give equivalent widths for 
comparison with recent work; agreement is within 
1-2 mA typically. Equivalent widths are missing 
where lines were not able to be measured due to 
"one-sided" blends or difficulty in defining the lo- 
cal continuum level. 

Notable in Table [T] is the discrepancy between 
Fe i and Fe n abundances. The abunda nces of 
Fe reported bv lAdelman fe Gulliver (|l990h do not 
show this dramatic lack of balance, and this might 
be viewed as supportin g the smaller temperature 
gradie nt d erived bv iGulliver. Hill, fe Adelman 



(|l994h and lHill. Gulliver, fe Adelmanl (|2004f ). In 
contrast with Fe, the abundances of Mg i and Mg 
ii shown in Table [1] do not show similar behavior. 
We discuss this result further below. 

5. Discussion 

5.1. How Does Rotation Affect Abun- 
dances? 




Fig. 4. — These plots show the derived abun- 
dances (data points) versus equivalent widths of 
Fe ii lines for two different assumed values of the 
microturbulence. Panels a) and b) show the de- 
rived abundances for microturbulence values of 
2 km s -1 and 4kms~ 1 , respectively. The dashed 
line shows the unweighted average abundances 
and the solid line shows the trend with equivalent 
width. We adopt a microturbulence of 2kms~ 1 in 
our abundance determinations. 



The main difference between a pole-on rapidly 
rotating star, as modeled here, and a classical 
plane-parallel stellar atmosphere model is that for 
the same integrated colors, the rotating model has 
some fraction of its surface at both higher and 
lower local effective temperatures than the non- 
rotating model. The expected effect is that there 
will be spectral lines in a range of excitation and 
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ionization energies where the two models give sim- 
ilar results. For Vega this is the case when the 
sum of the excitation energy and ionization en- 
ergy (for lines of ions such as Fe n) is about 10 eV. 
But for both higher and lower energy features the 
expanded range of temperatures will enhance line 
strengths, resulting in a decrease in the deduced 
abundances fro m those lines (described as "in- 
tensi fication" by Takeda, Kawanomoto. fe Ohishil 
20081) . 

We see that trend here. Table [2] summa- 
rizes the results of two r ecent analyses of Vega 
with non-rotating models dSadakane fe Nishimural 



ment, involves the ionization balance between Fe i 
and Fe n. Problems with the Fe i/Fe n ioniza- 
tion balance ha ve been reported for a wide range 
of stars (e. g., Gigasl 19861: lAllende Prieto et al 



19811 : lAdelman fe Gulliverl! 1990l ) along with the 
element-by-element results determined here. Our 
results from Mg i, Mg n, Al n, Si n, Ti n, Cr u, 
Mn i, and Fe n are about -0.6 dex below solar, 
roughly that found by earlier authors. These are 
mostly the dominant ionization stages and thus 
yield fairly stable abundances. Other ions, includ- 
ing C i, Ca i, Sc ii, Fe i, Ni i, and Ba n are deficient 
by about -1.0 dex, or even more, than the solar 
abundances (C i by -0.6), and are depressed by 
typically several tenths dex compared to the ear- 
lier studies, in this case because of their relatively 
small excitation and ionization energies. 

On the other side, the He abundance we de- 
duce, N Hc /N tot = 0.072 ± 0.004, which is essen- 
tially solar (0. 078), is substantially highe r than 
that found by lAdelman fe Gulliver (|l990h . run- 
ning counter to expectations. The result deter- 
mined here is based on five of the six lines in the 
ELODIE spectral window that are expected to be 
measurable (A 5875 is heavily involved with atmo- 
spheric water vapor lines), whi l e the helium abun- 
dance of lAdelman fe Gulliverl (l990) is based on 
A 4471 only. Otherwise, we have no explanation 
for why we obtain a larger helium abundance. 

In broadest terms, we find that if one can de- 
termine abundances from lines of the dominant 
ionization stage of an element, the errors induced 
by not accounting for rotation are small. Where 
lines from the dominant ionization stage are not 
accessible (e. g.,Ba in), one can expect large cor- 
rections to be required when standard, model- 
atmosphere analyses are applied to objects rotat- 
ing near breakup. 

One interesting example of the problems that 
can arise because of the corrections required be- 
tween different ionization stages of the same ele- 



19991 : iThevenin fe Idiardll999HJolrnsonll2002l ). For 

Vega, departures from LTE are predicted to pro- 
duce about 0.3 dex errors in abundances deduced 
from Fe i lin es while Fe n lines are barely affected 
(|Gigaslll986l ). However, these calculations are dif- 
ficult owing to the complexity of the atom and 
the lack of accurate collisio n and photoion i zation 
cross sections. For example, IPradhan et all (|l995h 
have found that many of the photoionization cross 
sections of Fe i are signific antly h i gher than those 
previously adopted (e. g., Gigasl 19861 ) with the 
possibility that the actual corrections from depar- 
tures from LTE are larger still. 

The problem with the Fe balance in Vega 
is confusing since at first glance, straightfor- 
ward LTE analyses dSadakane fe Nishimural 1981 



lAdelman fe Gulliverl 19901 ) provide apparent agree- 
ment between the abundances deduced from the 
two ions. This is in contrast to the sizable depar- 
tures from LTE required in other similar objects. 

However, even though we assume LTE in our 
analysis here, we also find a serious iron ioniza- 
tion imbalance amounting to ~ 0.4 dex, but in the 
opposite sense of that induced by non-LTE. To 
understand the origin of this imbalance we reana- 
lyzed representative lines from the two iron ioniza- 
tion states and, as a check, from the two magne- 
sium ions present, using a standard plane-parallel 
model. We find that rotation induces an apparent 
0.35 dex error in the Fe i/Fe n ionization balance, 
while the corresponding effect in the Mg i/Mg n 
balance is only about 0.1 dex. 

Thus we reach the amusing conclusion that a 
simple LTE analysis of Vega using models which 
do not account for rotation give a good ionization 
balance because of a nearly complete cancellation 
of the effects of photoionization-driven departures 
from LTE in the Fe i ion, on the one hand, and 
an enhanced Fe i line spectrum contributed by 
the extensive cool equatorial regions of the model 
owing to the favorable viewing geometry, on the 
other. Note however, the near balance between 
these two effects may disappear when one analy- 
ses lines in either the ultraviolet or infrared, owing 
to the changing relative contribution of the equa- 
torial regions to the overall light. 
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We note that Takeda, Kawanomoto. fc Ohishil 
(|2008l ) have independently commented on the near 
cancellation of departures from LTE versus the 
effects of rotation in the iron ionization balance. 
However, in their calculation the rotation induced 
errors are predicted to be about half those calcu- 
lated here, owing to the much lower rotation ve- 
locity and the corresponding dramatically reduced 
temperature gradient (~ 900 K) in their model. 

5.2. Is Vega a ABootis Star? 



Since Baschek fc Slettebak ( 19881 ) remarked 
that Vega showed an abundance pattern simi- 
lar to the ABootis stars, several stud ies (e. g., 
Venn fc Lambertlll99d : llliiic et al.|[l998h reported 
that Vega may be a mild ABootis star. We con- 
firm that result here. The abundance pattern we 
deduce matches well the main characteristics of 
the abundance patterns of ABootis st ars as sum- 
marized, for example by iHeiterl (|2002l ). Elements 
such as Si, S, Ca, and Sc fall in the middle of 
their respective typical ranges while O, Mg, Ti, 
Cr, Mn, and Fe are on the high side of normal 
and Ni and Ba are on the low side. While most 
elements fit the ABootis abundance pattern well, 
C and Al are somewhat out of the reported range. 
The Al abundance is based on one line and is 
not certain, while carbo n is off the lo wer end of 
the pattern r ei jorte d by IHeiterl (|2002h . However 
Pauzen et alj (|l999t ). in an extensive discussion of 
carbon and oxygen in this group of objects, find 
several objects with carbon abundances as low as 
-0.7 dex with respect to the Sun. We conclude 
that Vega would not be rejected as a A Bootis star 
on the basis of its carbon abundance and the rest 
of the abundances determined here are very much 
in keeping with membership in this group. 

5.3. Is Vega Well Mixed? 

A presumption, often unstated, about the na- 
ture of A Boo stars is that the deviations from 
sol ar composition are limited to surface layers (e . 
g.. iBaschekl Il992t iHolweger fc Sturenburgfll993f ). 
much the same as has been concluded for the Ap 
and Am stars which also occupy this part of the H- 
R diagram. But there has always been some con- 
cern about that assumption since unlike the latter 
groups the A Boo stars appear to have a distribu- 
tion o f rotation velocities similar to normal stars 
(e. g., Holweger fc Sturenbure 1993h . 



We argue here that since Vega is rotating at 
a significant fraction of breakup and yet displays 
fairly typical A Boo characteristics, it is unlikely 
that these composition anomalies are limited to 
the surface; more likely, Vega is well mixed. The 
literature on rotationally induced mixing has gen- 
erally focused on the surface layers and the ques- 
tion of whether the Ap and Am phenomena could 
be understood as due to diffusive separation (e. 
g., Charbonneau 19931 ) and not on how fast an in- 
homogeneity introduced on the surface would be 
mixed throughout the envelope. 

However, recent efforts to include the effects 
of rotation in evolutionary calcula tions of mas- 
sive stars (jMevnet fc Maederl Il997l) have led to 



an examination of how inho mogeneities will be 



redistributed throug h a star (jTalon et al.l 11997 



lEkstrom et aU 120081 ). suggesting that extensive 



mixing is to be expected. In fact, at the highest ve- 
locities in models down to 3 M Q , the lowest mass 
examined, the mixing is predicted to be so deep 
there is the possibility that some of the nuclear 
products from the CNO burning region might be 
mixed to the surface. 

This is an interesting possibility, given the low 
carbon abundance we have found. From this 
point of view, missing is an estimate of the ni- 
trogen abundance, the lines of which are out of 
the ELODIE spectral range. However measure- 
ments of nitrogen line equivalent widths have been 
reported elsewhere. To fill in the abundance 
of this important nuclide taking full account of 
the effects of rotation, we have calculated the 
abundanc es for the nitrogen equiva lent widths re- 
ported in IVenn fc Lambert! (jl990h for A 7442.28 
(log N/N tot = -4.05, [JV/JV to t] = +0.07) and 
A 7468.29 (log JV/JV tot = -4.02, [JV/JVtot] = +0.1) 
(A 7423 appears to be blended and we exclude it 
her e), finding values quite c lose to those deduced 
by IVenn fc Lambertl (|l990l ) at about 0.085 dex 
above solar. This is an intriguing result. Although 
it is difficult to know what "normal" is in this star, 
normalizing to oxygen gives [N/O] ~ +0.2 and 
[C/O] ~ -0.5, which may very well indicate that 
some CN cycle processed material has been mixed 
into the envelope of the star. In this regard we 
note that Vega represents a rather unique object; 
a few other A Boo objects have projected veloci- 
ties in the vicinity of 200 km s -1 , but Vega is the 
one object known to rotate as fast as 275kms~ 1 , 



10 



less than 10% from breakup in terms of angular ve- 
locity. However, without a better understanding 
of the composition of the material Vega started 
with, or other supporting information, we must 
leave this as just an intriguing possibility. 

In summary, we believe a fairly strong case 
can be made for the outer layers of Vega being 
well mixed, possibly even down to the edge of its 
nuclear burning core. If this is so then we are 
looking at about 2 M Q of material of highly un- 
usual composition in an object that is much too 
young to display such extreme depletion in heavy 
elements. In this case the various mechanisms 
put forth to explain the A Boo phenomena that 
rely o n its being limited to th e superficial layers 
(e. g.. iKamp fc Paunzenl [20021 ) seem excluded for 
Vega. Some form of dust - gas separa tio n, such as 
sugges ted bv lVenn fc Lambert! (|l990h or Hol weeerl 
may be involved but if so the mechanism 
likely must work at the time of Vega's formation 
since so much mass is involved. 

5.4. Determination of the Age and Mass 
of Vega 

We estimate the mass and age of Vega by lo- 
cating its measured luminosity and polar radius 
in an appropri a te evol ution grid, as described in 
Peterson et al ] lj2006br ). The interior models we 
adop t are from the BASTI databas^l (jPietrinferni 
20061 . and references therein) which include evolu- 
tionary calculations using scaled solar and alpha- 
enhanced compositions for stellar masses up to 
2AMq. The composition found here is not a 
perfect fit to either of the those mixtures, but 
the large enhancement of oxygen is about the 
same compared to the heavy metals as the alpha- 
enhanced mixture adopted there. Missing are the 
other alpha-rich elements at the enhanced levels, 
but given the dominance of oxygen even among 
these nuclei, that grid should give results more 
than adequate. To quantify how much the mis- 
match in the details of the distribution of abun- 
dances might affect the estimate we also calcu- 
late the mass and age using the scaled solar grid. 
In both cases, the heavy element fraction used is 
Z = 0.0093t^oo5 as calculated from Table Hand 
assuming [N e i/N to t] = —0.7 for abundances not 
obtained here. 



For the alpha-enhanced composition we obtain 
2.09 ± 0.03M© and 536 ± 29Myr for the mass 
and age. The simple scaled solar abundances in 
turn yield 2.14M© and 541 Myr and since the 
alpha-enhanced models are a much closer match 
to Vega's composition it is clear the errors intro- 
duced by the slight mismatch are small compared 
to the other uncertainties. 

As described in Peterson et al. I (|2006bh . whether 
Vega is solar composition throughout or the de- 
rived abundances represent the actual overall com- 
position, results in quite different estimates for 
the star's mass and age. Most previous authors 
have assumed an underlying solar composition 
yielding estimates of 2.3M© for the mass and 
an age in the neighborhood of 360 Myr. Since 
there is a distinct possibility that the composition 
we have derived applies to the star as a whole, 
Vega's estimated mass may be reduced and its 
implied age increased substantially. One immedi- 
ate consequence of this is a growing clash with the 
properties of the so-called " Castor moving group" 
( Barrado v Navascued 19981 ). which includes Cas- 
tor (a Gem), Fomalhaut (a PsA) and Alderamin 
(a Cep), in addition to Vega, and whose members 
are estimated to have an age of 200 ± 100 Myr. 
Even with an assumed solar composition Vega's 
age was not a comfortable fit for inclusion in this 
group. The increased age we propose would make 
it an unlikely member. 

This research was supported in part by a grant 
from the Naval Research Laboratory to D. M. P. 
and in part by NSF grant 06-07612 to Dr. Michal 
Simon. We also thank Dr. Robert L. Kurucz for 
extensive discussions. 
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Table 1 

The Abundance Analysis of Vega 



a 

(A) 



EP 

(cm - 



w\ logg/ log -jj^L 

mA) 



Blends 



Ho i ( Hc 

4471.498 

4713.139 

4921.931 

5015.678 

6678.154 

C i (log $ 

4770.021 

4771.730 

4775.889 

4932.050 



=0.072 ± 0.004 ) 



169087.008 




0.052 


0.070 




169086.864 


5 


-1.233 


0.078 


Fo n A 4713.193 


171135.000 


8 


-0.435 


0.060 




166277.546 




-0.820 


0.078 


Fo ii A 5015.755 


171135.000 


5 


0.329 


0.070 


Fc ii A 6677.306 


'— = -4.14 ±0.04, [JV c /JV t , 


3 ,] b = -0.62) 






60352.639 


7 


-2.052 


-4.16 




60393.148 


25 


-1.488 


-4.16 




60393.148 


7 


-2.013 


-4.16 




61981.818 


16 


-1.574 


-4.06 





5329.099 


■ = -3.32 ± 0.04, 


[N /N 


tot ] = -0.11) 




86625.757 


34 c 


-1.730 


-3.31 


5329.690 


86627.778 




-1.410 


-3.31 


5330.741 


86631.454 


24 


-1.120 


-3.31 


6046.438 


88631.146 


10 


-1.675 


-3.26 


6155.971 


86625.757 


77 d 


-1.051 


-3.36 


6156.778 


86627.778 




-0.731 


-3.36 


6158.187 


86631.454 


59 


-0.441 


-3.36 



Mg i 


(log 


^fe = -5.12 ± 

JVtot 


0.05, [JV Mg 


/Jv t ot] = -o.e 


,6) 




4702 


.991 


35051.264 


29 


-0.666 


-5 


.06 


5167 


.321 


21850.405 


81 


-1.030 


-5 


.06 


5172 


.684 


21870.464 


102 


-0.402 


-5 


.16 


5183 


.604 


21911.178 


119 


-0.180 


-5 


16 


5528 


.405 


35051.264 


28 


-0.620 


-5. 


.16 



Fe i A 5167.488 



Mg ii (log 



■'.'18 
«to 



-5.06 ± 0.04, [JV Mg /JVtot] = -0.6) 



4427.994 


80619.500 


-1.210 


-5 


.06 


4433.988 


80650.020 


-0.910 


-5 


.11 Fo i A 4433. 782 


4481.126 


71490.190 


0.740 


-5 


.01 



Al ii (log 
4663.046 



"Al 
JVtot 

85481.350 



6.22, [JV A i/JVtot] 



: -0.65) 
-0.284 



-6.22 



Si ii (log 


■jJ^L = -5.15 ± 0.05. 


, [JV SI 


/JVtot] = 


-0.66) 




4128.054 


'°79338.500 


32 


0.316 


-5 


.19 


4130.872 




79355.020 


54 


-0.824 


-5 


.19 


5055.984 




81251.320 


60 


0.593 


-5. 


.19 


6347.109 




65500.470 


118 


0.297 


-5 


.09 


6371.371 




65500.470 


82 


-0.003 


-5. 


.09 




= -5.01, [JVs/iVtot] = 


-0.3) 






6052.674 




63475.051 


7 


-0.740 


-5 


.01 


Ca I (log 


= -6.72 ± 0.12 


, [JVca/jVtot] = 


-1.04) 




4226.728 




0.000 




0.243 


-6 


.73 


4434.957 




15210.063 




-0.029 


-6 


.73 


4585.865 




20371.000 


1 


-0.386 


-6 


.68 


5588.749 




20371.000 


1 


0.210 


-6. 


.63 


5594.462 




20349.260 


5 


-0.050 


-6. 


.63 


5598.480 




20335.360 


4 


-0.220 


-6 


.63 


6162.173 




15315.943 


9 


0.100 


-6. 


.98 



Mn ii A 4128.129 



Mg ii A 6346.742 
Mg ii A 6346.964 



Fe i A 5598.287 



Sc ii (log 



'So 
Vtot 



-9.97 ± 0.05, [JVsc/jVtot] = -1.1) 
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Table 1 — Continued 



A 


EP 




logs/ 




Blends 


(A) 


(cm- 1 ) 


(mA) 






4246.822 


2540.950 


5 


0.320 


-10.02 




5526.79 


14261.320 


9 


0.130 


-9.92 




Ti n (log 


^EL = -7.65 ± 0.09, [Nt 


i/iV t ot] = 


-0.63) 




4468.507 


9118.260 


70 


-0.600 


-7.82 




4529.474 


12676.970 


9 


-1.830 


-7.69 




4563.761 


9850.900 


57 


-1.010 


-7.51 




4589.958 


9975.920 


16 


-1.790 


-7.61 


Cr ii A 4589.901 


4708.665 


9975.920 


3 


-2.410 


-7.69 




4779.985 


16515.860 


12 


-1.420 


-7.59 




4805.085 


16625.110 


21 


-1.100 


-7.59 




5336.771 


12758.110 


12 


-1.700 


-7.72 





Cr II (log 


= _ 6 . 91 ± o.l 


, [N C r 


Mot] = 


- 0.54) 


4252.632 


31117.390 


6 


-2.018 


-6.97 


4261.847 


25033.700 


18 


-3.004 


-6.92 


4554.988 


32836.680 


20 


-1.430 


-6.87 


4558.650 


32854.310 


61 


-0.660 


-6.87 


4565.740 


32603.400 


7 


-1.910 


-7.07 


4588.199 


32836.680 


48 


-0.830 


-6.87 


4592.049 


32854.950 


18 


-1.420 


-6.87 


4616.629 


32844.760 


16 


-1.530 


-6.87 


4618.803 


32854.950 


36 


-1.070 


-6.87 


4634.070 


32844.760 


29 


-1.220 


-6.82 


4812.337 


31168.580 


6 


-1.930 


-7.07 


4824.127 


31219.350 


39 


-1.220 


-6.72 


5334.869 


32844.760 


10 


-1.562 


-7.07 



Cr ii A 4261.913 



Mn i (log ■ 
4783.405 



"tot 

18531.663 



[N M n/N tot ] = -0.8) 
2 0.042 



-7.45 



Fo i (log - 



j^- = - 5.51 ± 0.1, [N Fo /N tot ] =-0.97) 
4132.058 °12968.553 29 -0.650 -5.54 



Fo i A 4131.935 
Fc i A 4131.971 



4134.677 


22838.321 


8 


-0.490 


-5.54 


Fc i A 4134.42 


4136.998 


27543.001 


4 


-0.540 


-5.54 




4250.119 


19912.494 


18 


-0.405 


-5.49 


Fc ii A 4250.437 


4250.787 


12560.933 


27 


-0.710 


-5.49 


Fc ii A 4250.437 


4260.474 


19350.890 


36 


-0.020 


-5.39 




4466.551 


22838.321 


9 


-0.590 


-5.36 




4476.019 


22946.814 


8 


-0.570 


-5.79 


Fc i A 4476.076 


4528.614 


17550.180 


15 


-1.072 


-5.51 




4918.994 


23110.937 


16 


-0.640 


-5.51 


Fc i A 4918.954 


4920.502 


22845.867 


27 


-3.955 


-5.51 


Cr II A 4920.23 


5324.179 


25899.987 


12 


-0.240 


-5.49 




5586.756 


27166.818 


11 


-0.210 


-5.59 


Fc ii A 5587.114 


5615.644 


26874.548 


15 


-0.140 


-5.44 




Fc ii (log 


-^Es. = .5.12 ± 0.09, [N Fl 


,/JVtot] = 


-0.58) 




4258.154 


21812.055 


14 


-0.467 


-5.34 


Fc ii A 4258.34 


4520.224 


22637.205 


45 


-2.990 


-5.13 




4522.634 


22939.358 


69 


-2.700 


-5.01 




4576.340 


22939.358 


23 


-3.390 


-5.06 




4582.835 


22939.358 


17 


-3.570 


-5.06 




4583.837 


22637.205 


88 


-2.490 


-4.97 


Fc ii A 4583.999 


4596.015 


50212.826 




-2.057 


-5.21 


Fc ii A 4595.682 


4620.521 


22810.357 


14 


-3.650 


-5.16 




4635.316 


48039.090 


13 


-1.650 


-5.21 
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Table 1 — Continued 



A 


EP 


w\ 


log gf 




Blends 


(A) 


(cm- 1 ) 


(mA) 







4656.981 


23317.633 


12 


-3.950 


-5.11 


Ti n A 4657.206 


4663.708 


23317.633 


5 


-4.145 


-5.11 




4666.758 


22810.357 


11 


-3.700 


-5.01 




4670.182 


20830.582 


8 


-4.350 


-5.11 


Sc ii A 4670.407 


4923.927 


23317.633 


114 


-1.820 


-5.06 




5534.847 


26170.181 


25 


-2.930 


-5.09 




6147.741 


31364.440 


14 


-2.721 


-5.24 




6149.258 


31368.450 


13 


-2.724 


-5.24 





Ni i (log 7^ = -6.79, [N m /N tot ] = -1.0) 

4714.417 °27260.894 3 0.160 -6.79 

Ba ii (log^ = -11.21, [JVna/ATtot] = -1-3) 

4554.029 ° 0.000 13 0.430 -11.21 

4934.076 0.000 7 -0.150 -11.21 FciA4934.005 



a For helium abundance, , r Hc 

"tot 

b [Jvel/J v tot ] = log^-log(^) Q 

c Thc equivalent width is for the blend with O I A 5629.690. 
d Thc equivalent width is for the blend with O I A 6156.778. 
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Table 2: Comparisons with the previous abun- 
dance studies 



Atomic 
Species 


log M/H 




[A r oi/A r tot] d 


SN a 


AG D 


This work 


He i 




-1.52 


-1.14 


-0.04 


C i 




-3.81 


-4.14 


-0.62 


N i c 






-4.53 


+0.09 


i 






-3.32 


-0.11 


Mg i 


-4.61 


-5.07 


-5.12 


-0.66 


Mg ii 


-4.96 


-5.11 


-5.06 


-0.60 


Al ii 




-6.33 


-6.22 


-0.65 


Si ii 






-5.15 


-0.66 


S i 






-5.01 


-0.30 


Ca i 


-6.11 


-6.21 


-6.72 


-1.04 


Sc ii 


-9.42 


-9.62 


-9.97 


-1.10 


Ti ii 


-7.31 


-7.47 


-7.65 


-0.63 


Cr ii 


-6.90 


-6.76 


-6.91 


-0.54 


Mn i 


-6.87 


-7.16 


-7.45 


-0.80 


Fe i 


-5.09 


-5.05 


-5.51 


-0.97 


Fe ii 


-5.09 


-5.12 


-5.12 


-0.61 


Ni i 


-5.94 


-6.38 


-6.79 


-1.00 


Ba ii 


-10.25 


-10.58 


-11.21 


-1.30 



Note. — The definition of the abundances we use differs from that adopted by SN and AG. For the helium abundance found 
here, the SN and AG abundances will be systematically larger than ours by 0.03 dex. 
" ISadakane fe Nishimural <198lfl 
lAdelman fe Gulliver! <|l990f l 
c Abundances from Table [T] 

d Solar abundances have been taken from lGrevesse fe Sauvall j!998f) 
e Abundance based on I Venn fe Lambert (|1990T > equivalent widths 
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